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We perform calculations relating the order parameter symmetry of organic quasi-one-dimensional 
superconductors to the bulk quasiparticle density of states and the bulk uniform spin susceptibility 
tensor at finite temperatures. Current experimental results suggest that some organic quasi-one- 
dimensional superconductors may exhibit triplet pairing symmetry. The purpose of our analysis is 
to attempt to narrow down the number of possibilities for the symmetry of the order parameter 
based on the current experimental evidence. 
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I. INTRODUCTION 

After the discovery of quasi-one-dimensional organic 
suuecconductorsEJ, many interesting aspects were stud- 
iedoS. The magnetic field versus temperature phase dia- 
gram for (TMTSF)2PF6 under pressure of 6 kbar was re- 
cently revisited by Lee et. aln. They found that the upper 
critical fields along the usual a, b', and c* directions were 
highly anisotropic. The upper critical field along a and 
b' displayed a strong positive curvature at lower temper- 
atures. Furthermore, the slope [~dHc2/dT]T^ for H || b' 



(a) 



> Hi^^ 

(b) 



was smaller than along the a direction, i.e., H, 
at high temperature. While, at low temperatures H, 
exceeded Hcf', after an unusual anisotropy inversion at 
H* = 1.6T. These initial results suggested the existence 
of triplet superconductivity in this system. An additional 
boost for the triplet scenario in (TMTSF)2PF6 ■was given 
by very recent NMR experimenta3'El. Lee et. alEB found 
that there is no ^^Se Knight shift in (TMTSF)2PF6 f( 
fields H II b' (P w 6 kbar)B, and H || a (P « 7 kbar)l 
These results suggest that Xb' ~ Xn, and Xa ~ Xn, 
where xn is the normal state susceptibility. Fjirthermore, 
a sharp and narrow (possibly Hebel-SlichterCI) peak was 
observed just below Tc in the ^"^Se NMR relaxatiiQu rate 



and P w 7 kbaiO. The 

suggests the existence of 



l/Ti, for H II a (iJ = L43, 
combined work of Lee et. al 
a triplet superconducting phase in (TMTSF)2PF6. How- 
ever, these recent ^"^Se NMR resuhs in (TMTSF)2PF6 
should be contrasted with earlier proton NMR results in 
(TMTSF)2C104eI which indicated the absence of the co- 
herence peak just below Tc — 1.03 K at = 0, and a 
behavior between Tc/2 and Tc. Based on their ex- 
periments, Takigawa et. aln argued that the supercon- 
ducting state of (TMTSF)2C104 has an anisotropic order 
parameter vanishing along lines on the Fermi surface, al- 
though they were not able to distinguish between singlet 
and triplet states. 

The temperature versus magnetic field phase diagram 
of (TMTSF)2C104|-at, ambient pressure was also mea- 
sured by Lee et. alB'i^, for H [| b'. Their results showed 
that the Pauli paramagnetic limit is also exceeded in 
this compound. Furthermore, Belin and Behnialld (BB) 



reported measurements of the thermal conductivity in 
the superconducting state of (TMTSF)2C104, indicat- 
ing that their data is inconsistent with the existence of 
gap nodes aX the Fermi surface as suggested by Taki- 
gawa et. aln. BB's afgument was that the behav- 
ior of the proton 1/TiB was limited only to T > Tc/2. 
Thus, it could not be considered convincing evidence for 
nodes in the gap because the temperature was not low 
enough. Even for conventionally gapped superconduc- 
tors, the exponential behavior of 1/Ti occurs only at very 
low temperatures (T <^ Tc). rJiietefore, these recent ex- 
perimental results combinedBtyO seem to suggest the 
existence of a fully gapped triplet superconducting state 
in (TMTSF)2C104. However, detailed ^^Se NMR exper- 
iments seem to be lacking for (TMTSF)2C104. 

Since both of these materials belong to the same class 
of compounds and have very similar crystal structures, 
it is feasible that they may exhibit similar order param- 
eter symmetries. However, additional experiments are 
necessary in order to make sure that the order param- 
eter symmetry of the (TMTSF)2C104 is the same as in 
(TMTSF)2PF6. A minimal set of experiments are nec- 
essary to identify the symmetry of the order parameter. 
For instance, data from upper critical field H || a, ^^Se 
Knight shift and NMR relaxation experiments are needed 
for CIO4, and information about the thermal conductiv- 
ity is needed for PFg. At zero magnetic field, however, 
these two materials have very different phase diagrams. 
The compound (TMTSF)2C104 is superconducting at 
low temperatures and zero pressure when slowly cooled 
(R-state), while (TMTSF)2PF6 is a spin density wave 
(SDW) insulator at low temperatures and low pressures, 
and becomes superconducting at low temperatures only 
at pressures above 5.5 kbar. Although the similar crystal 
structure of these systems suggests, from a simple group 
theoretical point of view, that the origin of the pair in- 
teraction is the same, the role of the proximity to an 
SDW phase in (TMTSF)2PF6 needs to be investigated 
both theoretically and experimentally. Furthermore, the 
true nature of the order parameter symmetry can only 
be directly evidenced via phase sensitiii&fixperimcnts like 
those performed in the cuprate oxideal3ll3. 

Abrikosov was the first to suggest the possibility of 
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triplet superconductivity in Bechgaard saltg0lll based 
on the behavior of (TMTSF)2C104 and (TMTSF)2PF6 
under X-ray bombardment. These systems exhibited a 
strong suppression of their critical temperature in the 
presence of radiatioiijnduced non-magnetic defectsE£lEZl. 
Gorkov and JerometS also suggested the possibility of 
triplet superconductivity based on a theoretical extrapo- 
lation of the semiclassical upper critical fields calculated 
near the zero field critical temperature. In addition, 
Lebedll3 pointed out that triplet superconductivity would 
manifest itself in Bechgaard salts through a remarkable 
reentrant phase in high magnetic fieMs. Later, Dupuis, 
Montambaux and Sa de Melo (DMS)Ea studied the field- 
versus-temperature phase diagram and the vortex lattice 
structure of these systems and concluded that these sys- 
tems could be either in a inhomogeneous singlet st ate U ke 
the Larkin-Ovchinikov-Fulde-Ferrel (LOFF) stateBE3 or 
in a triplet state. Furtiier studies of the LOFF state were 
performed by Dupuia23, while an equal spin triplet P^- 
ing state (ESTP) was discussed further by Sa de Melocj. 

After tha-iipper critical field measurements of Lee et. 
al.a, LebedEa was able to show that the LOFF state was 
Pauli paramagnetically limited, thus giving further sup- 
port for the triplet scenario (at least at high magnetic 
fields). The possibility of a magnetic field induced sin- 
glet (low fields) to triplet -(Jaich fields) transition was 
considered by Sa de MeloESEZl and Vaccarella and Sa 
de MeloEs, but currently there is no experimental ev- 
idence of a kink in the upper critical field of these sys- 
tems. Furthermore, reentrant superconductivity was^till 
expected for the triplet state at high magnetic fieldsEScEl. 
All these previous theories neglected the effects of fluc- 
tuations at high magnetic fields. A first attempt to in- 
corporate fluctuation effiects was made recently by Vac- 
carella and Sa de MeloE2l who showed that phase fluctua- 
tions can suppress the reentrant phase at high magnetic 
fields in the ESTE_state. Very recently Lebed, Machida 
and Ozaki (LM0)l3 suggested the possibility of a triplet 
phase with a d-vector with zero-component along the b' 
axis and a finite component along the a axis. A direct 
consequence of LMO's proposal is the existence of an 
anisotropic spin susceptibility at zero temperature and 
low magnetic fields: Xf = Xn, and Xa < Xw, where Xa 
corresponds to H || a, Xb' corresponds to H || b', and xn 
correspon(ia-|to the normal state susceptibility. However, 
Lee et. aluU reported new ^^Se NMR (Knight shift) re- 
sults for (TMTSF)2PF6 under pressure, indicating that 
Xb' « Xn, and Xa Xn- A fully gapped singlet "d- 
wave" order parameter for (TMTSF)2C104 was proposed 
by ShimaharaEil, while gapless triplet "f-wave" supercon- 
ductivity for (TMTSF)&PF6 was proposed by Kuroki, 
Arita, and Aoki_.(KAA)L3. Duncan, Vaccarella, and Sa 
de Melo (DVS)E3 performed a group theoretical analysis 
and suggested that a weak spin-orbit fully gapped triplet 
"Paj-wave" order parameter would be a good candidate 
for superconductivity in Bech gaard s alts. For this sym- 
metry Xa « XJV, ^'^fi^b' ~ XJvE^E^El in agreement with 
experimental resultsOij. In addition, this order parameter 



choice is consistent with the expectation of weak atomic 
spin-orbit effects, given that the heaviest element in these 
systems is Se. 

This paper is a longer version of the brief report by 
DVSlj. Here, we are concerned mostly with the sym- 
metry of the order parameter of a triplet quasi-one- 
dimensional superconductor at zero magnetic field, and 
have three main goals. First, we perform a group the- 
oretical analysis of the possible symmetries of the order 
parameter for an orthorhombic quasi-one-dimensional su- 
perconductor. Second, we calculate the temperature de- 
pendence of the bulk quasiparticle density of states and 
of the bulk uniform spin susceptibility tensor for vari- 
ous candidate symmetries of the order parameter consis- 
tent with the group theoretical analysis. Third, we make 
connections to scanning tunneling microscopy (STM) of 
quasiparticle density of states and Knight shift measure- 
ments of the spin susceptibility tensor. Thus, the re- 
mainder of the paper is organized as follows. In sec- 
tion II, we discuss the Hamiltonian used for weak spin- 
orbit (strong spin-orbit) coupling cases for generic sin- 
glet (pseudo-singlet) and triplet (pseudo-triplet) states. 
In section HI, we perform a group theoretical analysis of 



the symmetry of the order parameter, both in the weak 
and strong spin-orbit coupling limits. In section IV, we 
discuss the quasiparticle density of states for different 
symmetries of the order parameter. In addition, in sec- 
tion ^ we analyse the uniform spin susceptibility tensor 
Xmn for several order parameter symmetries, and we dis- 
cuss finite magnetic field effects and the role of vortices 
on Xmn ■ Finally, we summarize our results in section [v| . 



II. HAMILTONIAN 

In order to analyze the possible different symmetries 
of the order parameter we study quasi-one-dimensional 
systems with a single band in an orthorhombic lattice 
and allow for singlet or triplet pairing. We consider the 
following dispersion relation 



\tx\cos{k^a) - \ty\cos{kyb) - |t^| cos(A:^c), (1) 



where \tx\ ^ \ty\ ^ \tz\. Note that this, notation is 
slightly different from the standard notatioin. The corre- 
sponding translation is \tx:\ — > 2ta', \ty\ — > 2tb; \tz\ — > 2tc. 
Furthermore, a, b and c in our notation correspond to 
the unit cell lengths along the crystallographic directions 
a, b', and c* respectively. In the limit of weak interac- 
tions and low densities these quasi-one-dimensional sys- 
tems exhibit a well defined Fermi surface which is open, 
being formed of two separate sheets which intersect the 
Brillouin zone boundaries in the ky and k^ directions. We 
work with the Hamiltonian 



H — Hkin 



Hi, 



(2) 
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where the kinetic energy part is 



(3) 



where ek is the dispersion defined in Eq. ^, /i is the chem- 
ical potential, and 

= ^ E E ^"/37^(k, ^')blp{K <l)b,s{k', q) (4) 

kk'q a./J'yS 

is the interaction part with 



k+q/2,a^k+q/2,;3' 



(5) 



where the labels a, /3, 7 and S are spin indices and the 
labels k, k' and q represent linear momenta. 

In the case of weak spin-orbit coupling and triplet pair- 
ing (5 = 1), the model interaction tensor can be chosen 
to be 

Vc,p^s{k, k') = To.i3^sVr{-k, k')0r(k)(/.Hk'), (6) 

where Tafs^g = v^/j • with Vap = {i(J<Ty)ap- In ad- 

dition, the interaction Vr corresponds to the irreducible 
representation T with basis function 0r(k) representa- 
tive of the orthorhombic group. In the case of strong 
spin-orbit coupling and 5 = 1 the interaction 

V^p^siKk') = VriKk') [$r(k) ■ v,^] [<fr(k') • /2 

(7) 

where the interaction Vp corresponds to the irreducible 
representation T with basis function vector $r(k) rep- 
resentative of the orthorhombic group. In both weak 
and strong spin-orbit cpiipling, we use either the equa- 
tion of . niot ion methodEj or the functional integration 
methodtaEa in the zero center of mass momentum pair- 
ing approximation (which corresponds to the BCS limit 
in weak coupling) to obtain the anomalous Green's func- 
tion 



A„/3(k) 



(8) 



and the single particle Green's function 



ILUn 



(9) 



where = £k — ^J■ is the dispersion ek relative to the 
chemical potential ^, 



is the quasiparticle excitation energy, and 



At(k)A(k) /2, 



(10) 



(11) 



where the order parameter matrix A(k) has elements 
Aapik). The expressions for the anomalous (Eq. ||) 
and for the single particle (Eq. ^ Green's functions are 
valid only in the unitary case where A'l'(k)A(k) is diag- 
onal. In this initial work, we do not discuss in detail 
the non-unitary case for triplet superconductivity. Us- 
ing the single particle and anomalous Green's fuJtt : | ti o n s 
defined above and standard many body methodsEJ-L^Ej, 
we derive the order parameter equation in terms of the 
anomalous Green's function as 

A„/3(k) = -TJ2 ypo^ysiK k')F^s{k'), (12) 

,k' 

while the number equation, that fixes the chemical po- 
tential /X, can be written as 



N : 



d ~ 

G(k, iw„) — G"^(k, iuj„ 



(13) 



The matrix G(k, iujn) is the 2x2 block matrix 



(14) 



with matrix elements which are in turn 2x2 matrices in 
spin space. Performing the Matsubara sums and traces in 
Eq. (y) and Eq. (|l|) we obtain the more familiar forms 



A„,(k) = W(k,k')^^tanh j , 



(15) 



^-E-'' = E[(i + f)/(^'^) 



l-f-]il~f{Eu)) 



(16) 



r 



for the order parameter and number equations, re- spectively. These two equations must be solved self- 
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consistently away from the strict BCS limit in order 
to accommodate particle-hole asymmetries and strong 
coupling effects which tend to shift the chemical poten- 
tial substantially away from the Fermi energy. Quite 
generally these two equations are correct even in the 
strong mucling (or low density) regime provided that 



A. General Order Parameter 



The general form of the order parameter is 



A(k) = i (A,,(k)i + At,(k)d(k) • a) ay, 



(17) 



III. POSSIBLE SYMMETRIES FOR THE 
ORDER PARAMETER 

We will now discuss the symmetry of the order param- 
eter Ap^. For this purpose, we consider an orthorhombic 
crystaED with a conventional symmetry normal state, i.e., 
we assume that the normal state does not break the full 
lattice symmetry. By contrast, normal states that break 
full lattice symmetry are spatially inhomogeneous. Ex- 
amples of such spatially inhomogeneous states are stripe 
phases that occur in high-T^, superconductors, magnetic 
normal phases (ferro and antiferromagnetic), and normal 
states in a magnetic field (where the magnetic translation 
group needs to be incorporated). These states, however, 
will not be considered in the analysis performed in this 
section. 

In the case of weak spin-orbit coupling, the most gen- 
eral conventional symmetry group for the normal state is 
G„ = S0(3) X Gc X U(l) X T, where S0(3) is the group 
of rotations in spin space, Gc is the crystal space group, 
U(l) is the gauge group and T corresponds to time re- 
versal symmetry. Here, rotations of the spin and spa- 
tial degrees of freedom are independent. Since the su- 
perconducting state breaks U(l) symmetry, the possible 
order parameters correspond to different irreducible rep- 
resentations of Gs = SO (3) X Gc. However, in the case 
of strong spin-orbit coupling the normal state symmetry 



group is Gn ~ Gc''-' x U(l) x T, where now Gc'''' is iden- 
tical to the the space group Gc except when a rotation is 
involved. In this case, rotations of the spin and spatial 
degrees of freedom are not independent, and the possible 
order parameters correspond to different irreducible rep- 
resentations of Gs = Gc''"' . Although spin is not a good 
quantum number when spin-orbit coupling is important, 
we will still classify the possible states as singlet or triplet 
with respect to a pseudo-spin space provided that the su- 
perconducting state does not break inversion symmetry. 
The singlet/triplet pseujdo-spin classification was first in- 
troduced by AndersonLj in the context of heavy fermions. 
We will use this classification here whenever we refer 
to strong spin-orbit representations. Considering quasi- 
one-dimcaisional superconductors of orthorhombic crystal 
structura23, the relevant crystallographic point group ia 
D2h, which has only one dimensional representationstJ 
(See Table I . 



(J) 



which must transform according to the one dimensional 
representations of the orthorhombic point group D2h, un- 
der the assumption that the order parameter does not 
break the crystal translational symmetry, i.e., the order 
parameter is invariant under all primitive lattice trans- 
lations. The matrix 1 is the identity and the function 
Asi(k) is symmetric (even) under the transformation 
k — !■ — k. The three-dimensional vector d(k) is antisym- 
metric (odd) under k — > — k, while the function Atr(k) 
is symmetric (even) under the transformation k — k. 

We are mostly interested in triplet states which do 
not break time reversal symmetry, and our analysis is 
confined initially to zero magnetic field. However, for 
completeness, we shall mention briefly singlet and time- 
reversal-symmetry-breaking triplet states. In Table ||, 
a summary of the possible order parameter symmetries 
consistent with the D^h point group is presented for sin- 
glet systems assuming weak spin-orbit coupling, while 
in Table [H a summary of the possible order parameter 



symmetries for time-reversal-symmetry-breaking triplet 
states is shown. Tables EV and ^ summarize the group 
theoretical analysis performed for the order parameter 
matrix A(k) for time-reversal- invariant triplet supercon- 
ductors in the weak and strong spin-orbit coupling cases, 
respectively. All the tables mentioned above include the 
state nomenclature, the order parameter for the respec- 
tive singlet or triplet symmetries, and the loci of zeros 
of the corresponding quasiparticle excitation spectrum 
i?k, when /i = /i — min[ek] is positive. In Table [V, the 



vector (0,0,1) is indicated up to an arbitrary rotation in 
spm space. In Table 0, the numerical coefficients A, B 
and C define a specific model for the strong spin-orbit 
interaction appearing in Eq. (Q). It is very important to 
emphasize here that the basis functions A(k),y(k) and 
Z(k) transform like k^, ky and fc^, respectively, under the 
crystallographic point group operations, but they can not 
be chosen to be equal to , ky and kz . The reason being 
that the normal state dispersion ek intersects the bound- 
ary of the Brillouin zone along the ky and kz directions. 
Thus, it is necessary to take into account the periodicity 
of the order parameter matrix A(k) and of the order pa- 
rameter vector d(k) in reciprocal (momentum) space. In 
the case of orthorhombic crystals with dispersion ek de- 
fined in Eq. (|l|), where \tx\ ^ \ty\ 3> \tz\ a Fermi surface 
exists (for weak attractive interaction or high densities) 
and does intersect the Brillouin zones, therefore the min- 
imal basis set must be periodic and may be chosen to be 
A(k) = sin(fca;a), Y{k) — sin{kyb) and ^(k) — sm{kzc). 
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TABLE I: Character table for the D2h point group. 



Representation 


E 


C| 


C| 


Cf 


i 








Basis 


Ala 


1 


1 


1 


1 


1 


1 


1 


1 


1 


Big 


1 


1 


-1 


-1 


1 


1 


-1 


-1 


XY 


B2g 


1 


-1 


1 


-1 


1 


-1 


1 


-1 


XZ 


B-ig 


1 


-1 


-1 


1 


1 


-1 


-1 


1 


YZ 


Alu 


1 


1 


1 


1 


-1 


-1 


-1 


-1 


XYZ 


Blu 


1 


1 


-1 


-1 


-1 


-1 


1 


1 


Z 


B2u 


1 


-1 


1 


-1 


-1 


1 


-1 


1 


Y 


B:iu 


1 


-1 


-1 


1 


-1 


1 


1 


-1 


X 



TABLE IL Time reversal invariant singlet states in an orthorhombic crystal, assuming weak spin-orbit coupling. 



State Residual Group Order parameter Asi(k) = Q [fi > Q) = (/i < 0) 

^Aig SO3 X D2h X T 1 none none 

^Big SO3 X D2{Ci) X I xT XY lines none 

^B2g SO3 xD2{Cl)x I xT XZ lines none 

^Bsg SO3 xD2{Cl)x I XT YZ lines none 



B. Excitation Spectrum Characteristics 



For the weak spin-orbit coupling (Table IV) the only 
candidate for weaker attractive interaction (Jl> 0) is the 
state '^i?3«(a), where the quasiparticle excitation spec- 
trum i?k is fully gapped. In the case of strong spin-orbit 
coupling (Table ^ there are three candidates for weaker 
attractive interaction {fi > 0), i.e., the states Aiu,Biu 
and i?2ti, where the quasiparticle excitation spectrum 
may be fully gapped. When /i > 0, the state Aiu is fully 
gapped only for ^4 7^ and for any value of B and C; 
the state Biu is fully gapped only for B ^ and for any 
value of A and C; and the state B2U is fully gapped only 
for C ^ and for any value of A and B. 

We note in passing that in the case of strong attrac- 
tive interactions or small densities, where fl < 0, the 
excitation spectrum i?k is fully gappedcO for all of the 
states involving Tables |^ and for examplei^XJiere- 



fore, if Field Effect Transistor (FET) systemsBH can 
be constructed to change the density of carriers (n) in 
quasi-one-dimensional organic superconductors, it may 
be possible to study systematically the carrier density 
dependence of physical properties in these systems. The 
most interesting theoretical case corresponds to the sit- 
uation where i?k has zeros for jl > {n > Uc), and i^k 
has a full gap for fi < {n < ric). Here Uc is the critical 
density where fl — 0, where the characteristics of excita- 
tion energy iJk changes in a fundamental way. This effect 
here should-correspond to a Lifshitz transition for super- 
conductorgla where low temperature properties like spin 
susceptibility, specific heat, thermal conductivity, would 
change their qualitative behavior depending on n being 
greater, equal or less than Uc- The Lifshitz transition 
should also be accompanied by a topological phase tran- 



sition in the ground statJl^, where the zero temperature 
electronic compressibility diverges a,t n = Uc- However, 
if the excitation spectrum i5k is fully gapped (like in the 
case of the ■^i?3„(a) state) for // > (n > ric), E'k would 
remain fully gapped for /i < (n < ric), and only a 
crossover should occur. 

We also note in passing, that the states shown in Ta- 
ble [H are non-unitary since they break time reversal 
symmetry. In this case the excitation spectrum has two 
branches and is given by 



E,,±^^e, + \Atrmmdikw±\ciik)\], 

where q(k) — id*(k) x d(k). For p, > {n > ric), the 
regions in momentum space where i^k- = correspond 
to surfaces, while the regions where i?k+ — correspond 
to lines or the null set, as indicated in Table III. Thus, 



for /i > (n > ric), the excitation spectrum is gapless. 
However, for p, < {n < ric), there are no regions in 
momentum space where -Bk± — 0, which means that 
the superconductor becomes fully gapped. This dramatic 
change in the excitation spectrum (from gapless to fully 
gapped) of time reversal breaking (non-unitary) triplet 
states should also correspond to a Lifshitz transition at 
finite tempecatures and a topological transition in the 
ground statecJ. Although, there is currently no exper- 
imental evidence that quasi-one-dimensional supercon- 
ductor break time reversal symmetry (at zero magnetic 
field), it seems that in the case of Strontium Ruthenate 
time reversal symmetry breaking states are potentially 
good candidates for the order parameteiCJ. Strontium 
Ruthenate is very different from Bechgaard salts since it 
is non-organic and has essentially a tetragonal quasi-two- 
dimensional structure. If Field Effect Transistor (FET) 
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TABLE III: Time reversal breaking triplet states in an orthorhombic crystal, assuming weak spin-orbit coupling. 



State 


Residual Group 


Order parameter d(k) 


Sk± = (A > 0) 


£k± = (A < 0) 




D^{E) X Da X I{E) 


{l,i,0)XYZ 


surface (-), lines (+) 


none 


'Bi„(6) 


D^{E) X DiiCI) X I{E) 


(l,i,0)Z 


surface (-), lines (-I-) 


none 




D^{E) X D2(C|) X 1(E) 


(i,i,o)y 


surface (-), lines (-I-) 


none 




DUE) X £>2(Cf) X I{E) 


{l,i,0)X 


surface (-), none (+) 


none 


TABLE IV: Time reversal invariant triplet states in an orthorhombic 


crystal, assuming weak spin- 


-orbit coupling. 


State 


Residual Group 


Order parameter d(k) E^ = (/i > 0) 


= (/i < 0) 




I>oo(Coc) X X I{E) X T 


(0,0, i)a:f^ 


lines 


none 


^Biu{a) 


Doo(Coo) X D2(C|) X I{E) X T 


(0,0,1)^ 


lines 


none 


^B2u{a) 


D^Coo) X DiiC^) X I{E) X T 


(0,0, 1)F 


lines 


none 


''Bsuia) 


Doc(Coo)xD2(Cf)xJ(£)xr 


(0, 0, 1)X 


none 


none 


TABLE V: Time reversal invariant triplet states in an orthorhombic crystal, assuming strong spin-orbit coupling. 


State 


Residual Group Order parameter d(k) 


= (A > 0) 


Bk = (A < 0) 




D2 X I{E) X T 


{AX, BY, CZ) 


none, points, lines 


none 


Biu 


D2{Ci)xI{E)xT 


{AY, BX, CXYZ) 


none, lines 


none 


B2u 


D2(C|) X I{E) X T 


{AZ, BXYZ, CX) 


none, lines 


none 


Bsu 


Z?2(Cf) X I{E) X T 


{AXYZ, BZ, CY) 


points, lines 


none 



spectively. 

For the weak spin-orbit states is not only useful to 
know where the orbital components of the order param- 
eter change sign, but it is also important to visualize the 
d- vector at the Fermi surface. In the complete absence of 
the spin-orbit coupling the d-vectors can point in any di- 
rection with equal probability. However Bechgaard salts 
are expected to have weak spin-orbit coupling, leading to 
the weak pinning of the d- vector along a prefeijied direc- 
tion. Based on recent ^^5*6 NMR experimentaJD it seems 
that the d-vector in (TMTSF)2PF6 is not pointing along 
a or b' direction at least for magnetic fields beyond the 
1 T range. Therefore, in Fig. || we choose the weak spin- 
orbit coupling d-vectors to be pointing along the c direc- 
tion. We also illustrate in Fig. ^ the d-vector structure 
for the case of strong spin-orbit coupling. The d-vector 
for the states Aiu, -82^, and B^u (see Table ^ are 
shown in Figs. ||a, ||b, |c, and ||d, ioi A = B = 1 and 
C — 0, respectively. Notice the textured structure of the 
d-vectors at the Fermi surface. 

As discussed in the introduction the combined experi- 
mental evidence.|af the upper critical field measurements 
of Leo— et. aLoBIld, thermal conductivity measurements 
of BBHJ, and the Knigh.t_shift and NMR relaxation ex- 
periments of Lee et. alBu, suggest that (TMTSF)2C104 
and/or (TMTSF)2PF6 may be fully gapped triplet super- 
conductors. Thus, we consider next the temperature de- 
pendences of the quasiparticle density of states (QDOS) 
and of the uniform spin susceptibility tensor for a few po- 



systemsocS can be constructed to change the density 
of carriers {n) in Strontium Ruthenate, it may be also 
possible: (a) to study systematically the carrier density 
dependence of physical properties; (b) to search for a 
possible Lifshitz transition at finite temperatures, and 
(c) to investigate a possible topological quantum phase 
transition in the ground state. 



C. Three Dimensional Views of the Order 
Parameter 

A three dimensional view of the Fermi surface and 
order parameter sign for possible weak spin-orbit cou- 
pling singlet states compatible with the orthorhombic 
(1)2/1) group are shown shown in Fig. |^ for the case 
where jl > 0. In Fig. ||(a), (b), (c) and (d) the signs of 
the orbital component of the order parameter are shown 
for the ^Aig ("s-wave"), ^Big {"dxy-wave"), ^i?2g {''dxz- 
wave"), ^Big ("dj^^-wave"), respectively. The nodes are 
indicated by the dark thick lines separating the order 
parameter signs. In Fig. ^ a three dimensional view 
of the Fermi surface and order parameter sign for pos- 
sible weak-spin orbit coupling triplet states compatible 
with the orthorhombic {D2h) group are shown for the 
case where p. > 0. In Fig. ||(a), (b), (c) and (d) the 
signs of the orbital component of the order parameter 
are shown for the ^Ai„(a) ( "/a;yz-wave" ) , ^Biu{a) ("pz- 
wave"), ^B2u{a) ("py-wave"), ^Bi„(a) ( "pa;-wave" ) , re- 
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FIG. 1: Three dimensional views of the sign of the orbital component of the order parameter. Various weak spin-orbit coupling 
order parameter symmetries are illustrated (See Table (a) ^Aig ("s-wave"); (b) ^-Big ("da;j,-wave"); (c) ^-629 {"dxz-wave"); 
(d) ^Bsg ("dyz-wave"). The thick lines correspond to the zeros of the orbital parts of the order parameter as well as to the 
excitation energy node structure (_Ek = 0) on the Fermi surface (/i > 0). 



tential triplet candidate states, which are fully gapped: 

a) the weak spin-orbit coupling state '^B^uia) of Table IV; 

b) the strong spin-orbit couphng states Aiu of table |V| 
with A 0. Other fully gapped triplet states correspond 
to the Biu and B2u of table with B ^ 0, and C ^ 0, 
respectively. We will not discuss these states any further 
here. However, we will make comparisons of the cho- 
sen triplet candidates with the results of reference singlet 
states ^Aig ("s-wave") and ^Big ( "d^^-wave" ) . 



IV. QUASIPARTICLE DENSITY OF STATES 

The bulk quasiparticle density of states (QDOS) for 
these different symmetries can be obtained from the sin- 
gle particle Green's function as 

N{uj) = ^-TrY^ImGafi{k,ioJri ^oj + id), (18) 



where Ga/3(k, ?w„) is defined in Eq. (^. Features of the 
QDOS shown in Fig. ^ could be measured, in principle, 
during STM or photoemission experiments, depending 
on experimental resolution and material surface clean- 
ness. (In the case of photoemission, radiation damage 
due to X-rays, can also be an issue). Although, to our 



knowledge, these experiments have not yet been success- 
fully performed in (TMTSF)2C104 and (TMTSF)2PF6, 
our discussion can serve as qualitative guides for the ex- 
traction of gaps and symmetry dependent features of the 
experimental results when they become available. In par- 
ticular, STM or photoemission measured gaps could be 
compared with gaps measured thermodynamically (e.g., 
thermal conductivitylHI or specific heat). However, a 
note of caution is in place here. Due to the expected 
unconventionality of quasi-one-dimensional superconduc- 
tors, surface quasiparticle bound states could be present 
in these systems. The presence of these surface quasipar- 
ticle bound states was first noticed by Buchholtz and 
Zwickangl £BZ)c3 in isotropic p-wave superconductors 
and by HuEj in layered, d-wave superconductortw-. Re- 
cerjtljy, Sengupta et. alJiB followed the ideas of BZcJ and 
HueS and discussed the existence of these surface quasi- 
particle bound states in quasi-one-dimensional supercon- 
ductors. The QDOS that we discuss is a bulk property, 
thus the contribution of these surface states is not re- 
flected in our calculated QDOS. 

We show in Fig. | the bulk QDOS at T = K, and 
T = 1 K for the singlet states ^Aig ("s-wave") and ^Big 
( "dajy-wave" ) as references. While in Fig. ^ we present 
the QDOS at T = K and T = 1 K, with = 1.5 K, for 
the states '^i?3„(a) (weak spin-orbit coupling) and Aiu 
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FIG. 2: Three dimensional views of the sign of the orbital component of the d-vectors. Various weak spin-orbit coupling 
order parameter symmetries are illustrated (See Table ^v| ): (a) ^Ai„(a) ("/a;^z-wave"); (b) ^Bi„{a) ("p^-wave"); (c) '^B2u{a) 
("Py-wave"); (d) ^Biuia) ("pa;-wave"). The thick lines correspond to the zeros of the orbital parts of the d-vectors as well as 
to the excitation energy node structure (i5k = 0) on the Fermi surface {p. > 0). 



(strong spin-orbit coupling) for various values of the con- 
stants A, B, and C. These constants reflect the strength 
and the anisotropy of the pairing interaction Va/s-ysCk, k') 
defined in Eq. |^ for the strong spin-orbit coupling case. 
The symmetry dependent features of the QDOS are man- 
ifested through the magnitude of the order parameter 
vector d(k). For the ^-B3„(a) state the magnitude of 
d(k) is 

|d(k)| cx |sin(fc,a)|, (19) 
while for the state it takes the form |d(k) 

oc ^Ja^\ sin(fcxa)p + B^j sin(fcy6)|2 + C^l sin(fc2c)|2. 

(20) 

In Fig. ^(b), (c) and (d) we show only the case corre- 
sponding to C = 0, where A > A = B, and A <^ B, 
respectively. Notice that Fig. ^(d) is nearly identical to 
Fig. |^(a) given that |d(k)| is essentially the same in this 
case. Furthermore, notice that while the position of the 
peaks in Figs, ^(a), (b), (c) and (d) are essentially the 
same {cUp « ±5.3 K), the corresponding gap sizes are re- 
spectively ujg = 2.27 K; 0.29 K; 1.59 K; 2.27 K. In Fig. 
we show the frequency dependence near the gap edge at 
T = for the ^B3,j(a) (weak spin-orbit coupling) state 
indicated in Fig. ^(a). The characteristic shape exhibited 



here can also be found in the gapped states Ai^ shown 
in Fig. ^(b), (c), and (d). In addition, the frequency 
dependence near the gap edge in these states obeys a 
square-root power law. 

Although STM measures current (I) versus voltage (V) 
characteristics and not QDOS, the I-V (or dl/dV) char- 
acteristic are related to the QDOS. In a STM experi- 
ment one has to consider also the contribution of surface 
bound states, specially at zero voltage bias-ipdjiere a large 
contribution is expected in some casesEJ'E3E3. However, 
gap sizes, finite energy peaks and the general shape of 
the bulk QDOS should be, in principle, identifiable in an 
STM experiment at finite voltage bias. 

It is also desirable to have photoemission experi- 
ments probing QDOS in the superconducting state of 
(TMTSF)2C104 and (TMTSF)2PF6, however the sur- 
faces and the bulk of these materials seem to be very 
sensitive to X-rays. Thus, it is possible that photoemis- 
sion may be not able to probe the QDOS without in- 
troducing disorder at the surface and the bulk of these 
materials. Even if STM and photoemission experiments 
could be successfully performed, these experiments alone 
cannot uniquely determine the symmetry of the order pa- 
rameter in triplet quasi-one-dimensional superconductors 
since the QDOS depends only on the magnitude of d(k). 
Therefore, we now turn out attention to the calculation 




FIG. 3: Three dimensional views of d-vectors at the Fermi surface (/i > 0) for various weak spin-orbit couphng order parameter 
symmetries (See Table (a) ^Ai„(a) ('%j,z-wave"); (b) ^-Bi„(a) ("p^-wave"); (c) ^B2u{a) ("p^-wave"); (d) ^Biu(a) ("p^- 
wave"). For definiteness the d-vectors are chosen to point along the c direction. 



of the spin susceptibility tensor, which explicitly depends 
both on the magnitude and direction of d(k). 

V. SPIN SUSCEPTIBILITY 

We obtain a general form of the spin susceptibility ten- 
sor given by 

(21) 



where we used Einstein's notation that repeated Greek 
indices indicate summation, and the four-vector = 
(q, ii/). The tensor 

{Pmn)af3j5 = 9mgn{o'm)ap{o'n)fS: (22) 

contains the Pauli spin matrices and the scaled gyromag- 
netic factors gi = gi/2. The index £ takes into account 
the possibility of anisotropics in the gyromagnetic factors 
due to spin-orbit coupling. The other two tensors that 
appear in Eq. ^ are 



^Q/375(q, iuj)= P ^ ^ fI^{^ - q, w - iuj)Fi3s(k, iu), 



(23) 



(24) 



which contain the Green's functions described in Eqs. (|^) q ^ 0, we obtain an expression for the spin suscep- 
and (H), respectively. In the limit of w — > and tibility of a triplet superconductor (including lattice and 
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FIG. 4: Three dimensional views of d-vectors at the Fermi surface (Ji > 0) for various strong spin-orbit couphng order parameter 
symmetries (See Table (a) ylii,; (b) -Bii,; (c) i32u; (d) Biu, for A = B = 1, and C = 0. Notice that the direction of the 
d-vectors change along the Fermi surface, thus producing a textured structure. 
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FIG. 5: The frequency dependence of the QDOS N{lj) for singlet symmetries at T = K and T — 1 K. States ^Aig ("s-wave") 
and ^Big {"dxy-wa,ve") are shown respectively in (a) and (b). The parameters used are \tx\ ~ 5800 K, \ty \ = 1226 K, \tz \ = 48 K, 
and fi — —4101 K, with Tc — 1.5 K. N{uj) is in inverse temperature units (K~^) and lo is in temperature units (K). 



particle- hole asymmetry effects), 

Xriin,2 (k)] , (25) 

k 

where the k-dependent tensors have the forms 

(k) = g„rgnX\\0^)R^ ^m(k)c?n(k), (26) 



(k)- (k) (Sran - Re d:„(k)d„(k)) , (27) 



with (in(k) — rf„(k)/|d„(k)|. Here, the parallel compo- 
nent is 



Xl|(k) = -2Ai^ 



2 d f{E^) 



(28) 
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FIG. 6; Frequency dependence of the QDOS N{uj) for triplet states at T = K and T = 1 K. In (a) the weak spin-orbit 
coupling state ^B^u is shown. In (b) strong spin-orbit coupling state Aiu is shown for A = 0.20, B = 1.40, C = 0; in (c) for 
A = 1.00, B = 1.00, C = 0; in (d) for A = 1.41, B = 0.10, C = 0. The parameters used are \U\ = 5800 K, \ty\ = 1226 K and 
\tz\ = 48 K, and fi — —4101 K, with Tc = 1.5 K. N{uj) is in inverse temperature units (K~^) and uj is in temperature units 
(K). 



while the perpendicular component is 
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FIG. 7: A closer view of the T — curve in Fig. ]ma) show- 
ing the frequency dependence of the QDOS N{lu) near the 
gap edge for the ^Bsu state. This feature extends over an 
energy range of about 0.20 K (17 fieV). N{u}) is in inverse 
temperature units (K~^) and lo is in temperature units (K). 



The expressions derived in Eqs. (|2|) to (|2|) include some 
particle-hole symmetry effects, and arc valid at all tem- 
peratures, within the BCS limits. They dp-pot include, 
however, standard Fermi liquid correctionsEll, since these 
corrections are expected to be small for (TMTSF)2PF6 
at lower pressures (near the SDW phase). In this pres- 
sure regime this compound behaves like a good Fermi 
liquid. However, the same can not be said at higher 
pressures, where deviations from standard Fermi liquid 
behavior were reported^. The expression for Xmn (0,0) 
in Eq. (^) also allows for anisotropics in the scaled gy- 
romagnetic factors gi, where £ = m,n. For instance, 
these anisotropics could already exist in the normal state 
of the system. Thus, it is important to measure experi- 
mentally such anisotropics to determine the strength of 
the spin-orbit coupling in the normal state. For the D2h 
orthorhombic group Xmn is diagonal, but in the pres- 
ence of spin-orbit coupling the normal state Xmn can 
have, in principle, all diagonal elements different from 

each other i.e., xiV xiV xiT xiV ■ In the com- 
plete absence of spin-orbit coupling ge = 1 for all i, and 

(N) (N) (N) 

Xll =X22 -X33 ^XN- 

In Fig. @, the theoretical uniform Xmn is shown only for 
the ^B^u and Aiu states. We do not discuss other strong 
spin-orbit coupling states -Bi„ and -B2U, which can also 
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FIG. 8: Plot of the theoretical uniform spin susceptibility tensor components Xii (triangles); X22 (circles); X33 (squares) at 
low temperatures. The weak spin-orbit coupling state ^B^^ is shown in (a); strong spin-orbit coupling state Ai^ is shown in 
(b) for A = 0.20, B = 1.40, C = 0; in (c) for A = 1.00, B = 1.00, C = 0; in (d) for A = 1.41, B = 0.10, C = 0. The parameters 
used are \t^\ = 5800 K, \ty \ = 1226 K, \t^\ = 48 K, and = -4101 K, with = 1.5 K. 



be fully gapped because the atomic spin-orbit coupling 
in Bechgaard salts seems to be small (Se is the heavi- 
est element). The triangles correspond to xii circles to 
X22, and squares to X33. The susceptibilities Xmn are 
scaled by their normal state_values Xnm- It is known 
experimentally (Knight shift )El that the spin susceptibil- 
ity of (TMTSF)2PF6 for H || b' {H = 2.38 T) is very 

close to xi'^''- Knight shift experiments with H || aQ 
{H = 1.34 T) also indicate that the spin susceptibility 

is very close to Xa^^ ■ However if experiments indicate 
that the normal state susceptibilities obey the relation 

xi^^ ~ Xb^'' ^ Xc^^ ~ XN, then this is strongly sugges- 
tive that spin-orbit coupling effects are small, as expected 
from the fact that (TMTSF)2C104 and (TMTSF)2PF6 
have reasonably light-atoms. Furthermore, in recent ex- 
periments, Lee et. aln observed a coherence peak in the 
NMR relaxation rate 1/Ti for H \\ a, H = 1.43 T. In 
conventional superconductors the enhancement of 1/Ti 
near Tc is often associated with a divergence, of the den- 
sity of states at the edge of the energy gapQ. Although 
the weak spin-orbit couphng fully gapped state ■^B^uia) 
( "pa;-wave" ) (see Fig. ^) does not exhibit a divergence in 
QDOS at the edge of the energy gap, the QDOS changes 
very rapidly at the edge and may produce a peak or a 
kink in l/Ti for ^^Se. 



These results reinforce the evidence for triplet super- 
conductivity in (TMTSF)2PF6. For fields along the c* 
direction the superconducting state is easily destroyed, 
which makes impractical the Knight shift measurement 
for ""^Se. This measurement requires magnetic fields of 
the order of one Tesla, which is an order of magnitude 
higher than the upper critical field along this directionQ. 
For definiteness, we choose the unit vectors 1 (m = 1), 2 
(to = 2) and 3 (m = 3), to point along the a, b', and c* 
directions, respectively. 



A. Finite Magnetic Field Effects 

In this section, we discuss qualitatively some of the 
effects of a finite magnetic field upon the spin suscepti- 
bility tensor. We consider first the effects of vortices on 
the spin susceptibility tensor within a simple two fluid 
model, and second the possible rotation of the d-vector 
in the presence of a finite magnetic field. 



1. Role of Vortices 

Here, we would like to comment on the role of vor- 
tices when the spin susceptibility tensor is measured. 
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The Knight shift measurements of Lee et. alBn in 
(TMTSF)2PF6 had to be performed at reasonably high 
fields H = 2.38 T for H || b and = 1.43 T for H || a in 
order to obtain a clear NMR signal. In these field regimes 
the superconductor (TMTSF)2PF6 is in an inhomoge- 
neous— vortex state. This vortex state jSKas analysed by 
DMScj, and later shown by Sa de MeloEil to correspond 
to a rectangular Abrikosov vortex lattice (in the semi- 
classical low field regime) and to a rectangular Josephson 
vortex lattice (in the quantum high field regime) for equal 
spin triplet pairing or singlet pairing. We will use this 
result to construct a simple phenomenological two-fluid 
model to describe the effect of vortices on the spin sus- 
ceptibility tepspr, while a more sophisticated calculation 
is under wayua. Let us consider our superconductor at 
finite T and H, and assume a semiclassical description 
of the vortex state. Assume that the magnetic field is 
applied along the h direction {Hn) and produces a rect- 
angular Abrikosov lattice with lattice constants and 
ip, where a and P are direction perpendicular to h. The 
density of vortices is Hy — 1/^a^^- Assuming that the 
vortices have finite core size, characterizedphy lengths 
and £,0, and that the vortex core is normalEHI, the normal 
fraction of the superconductor is 



(30) 



Vortex core states have -b|ipji^-jCp|tmi| Siyc} i^ investigated 
in singlet superconductord^Hl^'EjOEilLJH'LJ; however, 
these states, in triplet superconductors remain to be in- 
vestigatedE3. We treat core states very crudely as nearly 
particle-in-a-well-Vike with characteristic sizes controlled 
by the coherence lengths ^a, and transverse to the 
magnetic field direction. This is a crude improvement 
upon treating vortex cores as purely normal, while a more 
sophisticated calculation of vortex core states is under 
wayc3. In the present case, if the two fluids are viewed as 
nearly independentEJ, then the total spin susceptibility 
tensor of the superconductor is the sum of the response 
of its individual components. Thus, we write 

XraniT,Hn) = /^^{T, Hn)nf + x'HiT, Hn)il - n/), 

(31) 



(C) 

where Xmn(T, Hn) is the vortex core contribution, and 
XmniT, Hn) IS the superconducting contribution. How- 
ever it is best to illustrate the qualitative behavior of 
Xmn{T^^ Hfi) via the ratio 

R{T, Hn) = Rc{T, Hn)nf + RgiT, Hn){l ~ nj) (32) 
Here R(T,Hn) = iJft)/xS(T, i/^, and 

Rj{T,Hn) = XmUr,i/ft)/xS(T,ifft), with j = C,S. 

In Fig. ^a the ratio R{T, Hn) for a weak spin-orbit cou- 
pling triplet superconductor is illustrated for fixed tem- 
perature T < Tc- It is assumed that the applied mag- 
netic field is perpendicular to the d-vector (H _L d). 



and that d-vector direction is not affected by the mag- 
netic field. In this case, < Rc{T,Hn) < 1, such that 
Rc{T,Hn) changes from a smaller value at small n/ at 
fixed temperature, and as n/ tends to 1, Rc{T, Hn) con- 
verges to 1, given that the core spectrum becomes the 
normal state continuum at Hn = Hc^n- The supercon- 
ducting term is assumed to be essentially constant for 
(H _L d), and Rs{T, Hn) « 1, implying that the total ra- 
tio R{T,Hn) « 1. These results are largely independent 
of temperature and field within the approximation used. 
However, the discussion above serves only as a qualita- 
tive|-ttuide, while a more sophisticated analysis is under 
wa5ica. In Fig. ||b the qualitatively similar situations of a 
weak spin-orbit coupling singlet superconductor or weak 
spin-orbit coupling triplet superconductor with H || d 
are illustrated. In this case, the core term Rc{T, Ha) still 
satisfies the condition < Rc{T,Hn) < 1, and grows at 
fixed temperature from a smaller value at small Uf to 
Rc{T, Hn) = I as rif tends to 1. In addition, the contri- 
bution of the superconducting part Rs{T, Hn) depends 
strongly on magnetic field (at fixed temperature) grow- 
ing from a small value at small to Rs(T, Hn) = 1 as 
— > 1. The condition < Rs{T,Hn) < 1 is also satis- 
fied. Therefore, R{T,Hn) would be very strongly depen- 
dent on magnetic field (at fixed temperature) in (a) the 
weak spin-orbit triplet case for H || d , or in (b) the weak 
spin-orbit singlet case. Thus, if measurements of tem- 
perature and magnetic field dependence of Xmn{T, Hn) 
in (TMTSF)2PF6 are such that Rs{T,Hn) « 1 (being 
largely independent of magnetic field at fixed tempera- 
ture), then a weak spin-orbit triplet state would be con- 
sistent with such finding for H ± d. It is possible, how- 
ever, that a small magnetic field can rotate the d-vector 
in the weak spin-orbit coupling limit, and this effect is 
discussed next. 



2. Rotation of d-vector 

In what follows, we will discuss the possible rotation of 
d(k), not including vortex core contributions, i.e., we will 

consider effectively only Xmn{T, Hn). However, similar 
effects are expected for the core contribution. 

For the orthorhombic symmetry Xmn is diagonal and 
is calculated under the assumption that d(k) is constant, 
i.e., the direction of d(k) is assumed not to change upon 
the application of a small magnetic field. In the case of 
state Aiu (strong spin-orbit coupling) a small magnetic 
field cannot rotate d(k) which is pinned to a particu- 
lar lattice direction. In the strong spin-orbit coupling 
case Xmn is still diagonal, however the diagonal compo- 
nents are not equal in general. Thus, the experimentally 
measured Xmn ^^'^ the theoretically calculated Xmn (^t 
constant d(k)) should agree for small enough magnetic 
fields (see Fig. ||b and ||d). However, in the case 
of state ■^i?3„(a) (weak spin-orbit coupling) a small mag- 
netic field can easily rotate d(k) to be perpendicular to H 
(barring any additional orbital effects due to the coupling 
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FIG. 9: Schematic plots of the ratios R (sohd hne); Rc 
(dashed line); and Rs (dotted line) are shown in (a) for the 
case of a weak spin-orbit coupling triplet superconductor with 
H _L d; and in (b) for the case of a weak spin-orbit coupling 
singlet superconductor or weak spin-orbit coupling triplet su- 
perconductor with H II d. The ratio R is connected to Rc 
and Rs via Eq. ^ (See text following Eq. ^ for definition 
of ratios R, Rc and Rs-) 



of H with charge degrees of freedom), thus minimizing 
the magnetic free-energy Fmag = -HrnXmnHn/2. In the 



weak spin-orbit coupling case, Xmn ~ XN^mm where xn 
is the normal state value, for any direction of H. Thus, 
what experimentalists measure, Xmm ^'^d what we cal- 
culate, Xmrai (shown in Fig. ||a) jaapuld be different. (This 
point was considered by LeggettcZl in his discussion of the 
phases of liquid ^He). Another important distinguish- 
ing feature between strong spin-orbit and weak spin-orbit 
coupling scenarios is the magnitude of the flop field Hf 
required to rotate the d-vector from a preferred lattice 
direction. In the case of strong spin-orbit coupling the 
flop fields are in the Tesla range, while for weak spin-orbit 
coupling the flop fields are in the Gauss range. For in- 
stance, based on the discussion above, the state ^B^u{a) 
(weak spin-orbit coupling "pa;-wave") would present no 
Knight shift for any direction of the externally applied 
magnetic field provided that the magnitude of H is larger 
than the small spin-orbit pinning field Hf for any given 
direction. Under these conditions, no Knight shift should 
be observed for fields along the a, b' or c*-axis or for any 
set of angles with respect to these axis. Furthermore, in 
the weak spin-orbit coupling state ^-63^(0) ("p^-wave"), 
the normal state susceptibilities should be nearly iden- 
tical, i.e. xii^^ ~ X22^^ ~ Xas'^ = Xn, since the scaled 
gyromagnetic factors ~ 1 for all 



FIG. 10: Scenarios for the behavior of the experimen- 
tal Xnn/Xnn Considering that the d-vector is pinned along 
the n direction. The ratio Xnn/x^S is schematically plot- 
ted as a function of applied magnetic field at a constant 
T = (left set) and T > (right set). The lower up- 
per critical field J/c""* is assumed to be very small and it 
is not shown. In each of these graphs, the left and right 
vertical dotted lines represent the minimum field {H^}^) re- 
quired to obtain Knight shift experimental data and the up- 
per critical field {H^"^) corresponding a fixed temperature 
T. Three possible fiop transition scenarios at T = 0: (a) 

min 



Hf 



< H'^l^ < nf (flop transition beyond present ex- 
perimental reach); (b) H^}^ < Hf < Hf (flop transition 
can be observed if jump in susceptibility is large enough); 
(c) Hf^ < Hf < Hf (fiop transition may be difficult to 
be observed since jump in susceptibility may be too small). 
Similar scenarios are illustrated in (d), (e) and (f) for T > 0. 



The spin susceptibility tensor for a triplet supercon- 
ductor is very sensitive to magnetic fields both in the 
weak or strong spin-orbit coupling scenarios since the 
applied field can depin the d-vector from a preferred lat- 
tice direction and cause a flop transition with a corre- 
sponding change in the spin susceptibility. This means 
that with increasing magnetic field at fixed temperature 
T, H > Hc^, and with H || d initially, the d-vector 
flops a± H — Hf < Hc2- Therefore, beyond Hf the d- 
vector becomes perpendicular to the applied magnetic 
field (d ± H). This flop transition produces a discon- 
tinuous change in the spin susceptibility from smaller to 
larger (normal) values, as the d-vector rearranges itself 
from d II H to d ± H. 

Consider for definiteness that the d-vector is pinned 
along a given lattice direction n, and that H || d. The 
observability of this flop transition depends on how the 
magnitude Hf^ compares to the magnitudes of the min- 
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imum field required to perform a Knight shift experi- 
ment, H^^^, the lower Hc^^ and the upper i/i""* crit- 
ical fields along the given direction n. Three possible 
scenarios are illustrated in Fig. ^ where Tji"-* is con- 
sidered to be very small and it is not indicated in the 
figure. If ^fj"'' "C H^^l^ then the transition cannot be 

observed in a Knight shift experiment. If -ffj"' ^ Hc2 ^ 
then, depending on experimental resolution, the transi- 
tion may not be resolvable and the data would seem char- 
acteristic of a singlet response. However, the flop transi- 
tion could be easily measured in the intermediate regime 
^min < ^/"^ ^ Hc2 ^ . If d was pinned along n — a, then 

for H II a such that H > 11^'' a flop transition could take 
place. Similarly situations would occur if the d-vectors 
were pinned along fi — h' (c*), then for H || b' (c*) such 

that H > h'^^ ^ [H > H^^) a flop transition could take 
place. So far there is no experimental evidence of this 
possible d-vector flop transition. 



3. Connection to Anisotropy Inversion in 

It is also important to make ja-connection between 
spin susceptibility jneasurementaju and upper crjijbical 
fleld measurementfjd in (TMTSF)2PF6. Lee et. alE ob- 
served an anisotropy inversion between -ffi^^ and-jffi^ 
at H* « 1.6 T. It was suggested theoreticallyE3 that 
this anisotropy inversion was related to the presence of 
a component of the d-vector along the a direction, such 
that -ffi^'' (H II a) would be paEamagnetically limited. In 
addition, it was suggested thatEj the d-vector would have 
zero component along the b' direction, such that Hc^ 
(H II b') would not be paramagnetically limited. While 
this suggestion is somewhat appealing, it seems to indi- 
cate that the d-vector has a strongly pinned component 
along the a axis. 

A comparison between the upper critical field experi- 
ments of Lee et. alB for (TMTSF)2PF6, witk^the upper 
critical field experiments of Shivaram et. a/.EJ for UPta 
can be illuminating. In UPta an anisotropy inversion oc- 
curs between the upper critical fields between iji^"* and 
H^;' at H* « 1.9 T. This anisotropy inversiqn in UPtg 
was successfully explained by Choi and SaulsH under the 
assumption that is paramagnetically limited, while 
Tfc^^ is not. The nearly perfect fitting of the experi- 
mental curves of UPts required the assumption that the 
d-vector was locked to the c axis by a strong spin-orbit 
coupling. Considering that the values of the anisotropy 
inversion field H* for (TMTSF)2PF6 and UPta arCjMjry 
similar in magnitude, the theoretical interpretationE3 of 
the anisotropy inversion in (TMTSF)2PF6 seems to as- 
sume implicitly a strong spin-orbit coupling, which is 
very unlikely for Bechgaard salts given that the heaviest 
element is Se. The ratio between (atomic) spin-orbit 
couphngs in (TMTSF)2PF6 (for Se) and UPtg (for U) 



can be estimated to be small: r^o ^ 0.15. This implies 
that the physical origin of the anisotropy inversion in 
(TMTSF)2PF6 should not be implicitly or exphcitly the 
same as in UPts. Furthermore, the theoretical sugges- 
tion connecting the anisotropy inversion with spin sus- 
ceptibility measurements^ would lead to Xa ~ 0-'^X% at 
T — Q. However, experiments indicate that Xa ~ Xnu in 
(TMTSF)2PF6 (Tc = 1.2 K) down to T = 0.32 K. 

It is important to note that the related compound 
(TMTSF)2AsF6 has an SDW phase below = 12 K 
with easy, intermediate andJiard axis along the b', a, and 
c* directions, respectively^. This system also presents a 
spin-fiop transition where the spin orientation fiops from 
the b' axis to an orientation predominantly parallel to 
sB. Given that the SDW phase in (TMTSF)2ES'6 is 
very similar to the SDW phase in (TMTSF)2AsF6B, we 
infer that the connection betwega. the experimental re- 
sults Xa ~ Xn and Xfc' ~ Xa^™ with the anisotropy 
inversion of the upper critical fields Hc2^ and Hc^ in 
(TMTSF)2PF6 may require a deeper understanding of 
the interplay between the spin-density-wave and the su- 
perconducting phases. Thus, it is important to search 
for an Hc2 anisotropy inversion in the sister compound 
(TMTSF)2C104 in order to have a more complete pic- 
ture of the superconducting state in the Bechgaard salt 
family (TMTSF)2X. To our knowledge detailed experi- 
mental studies of Hc2^ (high fields and low temperatures) 
have not yet been performed for (TMTSF)2C104. Hav- 
ing made this last experimental connection, we are ready 
to summarize our results. 



VI. SUMMARY 

We performed a group theoretical analysis of the possi- 
ble symmetries compatible with the D2h (orthorhombic) 
point group, and focused on the weak and strong spin- 
orbit coupling triplet cases at zero magnetic field. We 
also discussed order parameter symmetry features and 
temperature dependence of the quasiparticle density of 
states and spin susceptibility tensor of ap, orthorhom- 
bic quasi-one-dimensional superconductorLll. Based on 
current experimental evidence and the assumption that 
the origin of superconductivity in (TMTSF)2C104 and 
(TMTSF)2PF6 is essentially the same, we suggested that 
the weak spin-orbit coupling state ^B^uia) ("Pa;-wave") 
is a very good candidate for the order parameter sym- 
metry for these systems since this state is: (1) fully 
gapped and. consistent with thermal conductivity mea- 
surementsEll; (2) characterized by weak spin-orbit cou- 
pling and consistent with weak spin-orbit scattering fits 
oiTc{H) for (TMTSF)2PF6N at low magnetic fields; (SA 
consistent with no observable Knight shift when H || aS 
or II b'cl, and predicted to have no observable Knight shift 
for any direction n of H provided that the magnitude of 
the applied field is larger than the small spin-orbit pin- 
ing field Hf\ We would like to thank NSF (Grant No. 



16 



DMR-9803111) and its REU program for support. One ^ Present address: Department of Physics, Mas- 
of us (C. A. R. Sa de Melo) would like to thank the Aspen sachusetts Institute of Technology, Cambridge, MA 
Center for Physics for its hospitality. 02139-4307. 



^ D. Jerome, A. Mazaud, M. Ribault, and K. Bechgaard, J. 
Phys. Lett. 41, L95 (1980). 

^ T. Ishiguro and K. Yamaji, Organic Superconductors 
(Springer- Verlag, 1989). 

^ J. M. Williams, J. R. Ferraro, R. J. Thorn, K. G. Carlson, 
U. Geiser, H. H. Wang, A. M. Kini, and M. H. Whang- 
boo, Organic Superconductors (Including Fullerenes) : Syn- 
thesis, Structure, Properties and Theory (Prentice Hall, 
1992). 

" I. J. Lee, M. J. Naughton, G. M. Banner, and P. M. 

Chaikin, Phys. Rev. Lett. 78, 3555 (1997). 
^ I. J. Loo, D. S. Chow, W. G. Clark, M. J. Strouse, M. J. 
Naughton, P. M. Chaikin, and S. E. Brown (2000), cond- 
mat/0001332. 

^ I. J. Lee, S. E. Brown, W. G. Clark, M. J. Strouse, M. J. 
Naughton, W. Kang, and P. M. Chaikin, Phys. Rev. Lett. 
88, 017004 (2002). 

L. C. Hebel and C. P. Shchter, Phys. Rev. 113, 1504 

(1959). 

® M. Takigawa, H. Yasuoka, and G. Saito, Journ. Phys. Soc. 

Japan 56, 873 (1987). 
^ I. J. Lee, A. P. Hope, M. J. Leone, and M. J. Naughton, 

Apphed Superconductivity 2, 753 (1994). 
^° 1. J. Lee, A. P. Hope, M. J. Leone, and M. J. Naughton, 

Synth. Metals 70, 747 (1995). 
" S. Belin and K. Bchnia, Phys. Rev. Lett. 79, 2125 (1997). 
12 D. A. Wolhnan, D. J. V. Harlingen, W. C. Lee, D. M. 

Ginsberg, and A. J. Leggett, Phys. Rev. Lett. 71, 2134 

(1993). 

1^ A. Mathai, Y. Gim, R. C. Black, A. Amar, and F. C. 

Wellstood, Phys. Rev. Lett. 74, 4523 (1995). 
" A. A. Abrikosov, JETP Lett. 37, 503 (1983). 
15 A. A. Abrikosov, J. Low Temp. Phys. 53, 359 (1983). 
1^ M. Y. Choi, P. M. Chaikin, S. Z. Huang, P. Haen, E. M. 

Engler, and R. L. Greene, Phys. Rev. B 25, 6208 (1982). 

S. Bouffard, M. Ribault, R. Brussetti, D. Jerome, and 

K. Bechgaard, J. Phys. C: Solid State Phys. 15, 2951 

(1982). 

1* L. P. Gorkov and D. Jerome, J. Physique Lett. 46, L643 
(1985). 

1^ A. G. Lebed, Sov. Phys. JETP Lett. 44, 89 (1986). 
2° N. Dupuis, G. Montambaux, and C. A. R. Sa de Melo, 
Phys. Rev. Lett. 70, 2613 (1993). 

21 A. 1. Larkin and Y. N. Ovchinikov, Sov. Phys. JETP 20, 
762 (1965). 

22 P. Fulde and R. A. Ferrel, Phys. Rev. 135, A550 (1964). 

23 N. Dupuis, Phys. Rev. B 51, 9074 (1995). 

2* C. A. R. Sa de Melo, Physica C 260, 224 (1996). 

25 A. G. Lebed, Phys. Rev. B 59, R721 (1999). 

2® C. A. R. Sa de Melo, The Superconducting State in Mag- 
netic Fields: Special Topics and New Trends (World Sci- 
entific, Singapore, 1998), chap. 15, pp. 296-324, (Edited 
by C. A. R. Sa de Melo). 

2^ C. A. R. Sa de Melo, J. Suporcond. 12, 459 (1999). 

2« C. D. Vaccarella and C. A. R. Sa de Melo, Physica C 341- 
348, 293 (2000). 



2® C. D. Vaccarella and C. A. R. Sa de Melo, Phys. Rev. B 

63, R180505 (2001). 
3° A. G. Lebed, K. Machida, and M. Ozaki, Phys. Rev. B 62, 

R795 (2000). 

31 H. Shimahara, Phys. Rev. B 61, R14936 (2000). 

32 K. Kuroki, R. Arita, and H. Aoki, Phys. Rev. B 63, 094509 

(2001). 

33 R. D. Duncan, C. D. Vaccarella, and C. A. R. Sa de Melo, 
Phys. Rev. B 64, 172503 (2001). 

3^* V. P. Mineev and K. V. Samokhin, Introduction to Uncon- 
ventional Superconductivity (Gordon and Breach, Amster- 
dam, 1999). 

35 C. A. R. Sa de Melo, M. Randeria, and J. R. Engelbrecht, 

Phys. Rev. Lett. 71, 3202 (1993). 
3® J. R. Engelbrecht, M. Randeria, and C. A. R. Sa de Melo, 

Phys. Rev. B 55, 15153 (1997). 
3'' The Bechgaard salts (TMTSF)2C104 and (TMTSF)2PF6 

have truly a triclinic lattice structure. 
3® P. W. Anderson, Phys. Rev. B 30, 4000 (1984). 
3® M. Tinkham, Group Theory and Quantum Mechanics 

(McGraw-Hill, New York, 1964). 

R. D. Duncan and C. A. R. Sa de Melo, Phys. Rev. B. 62, 
9675 (2000). 

*i J. H. Schon, C. Kloc, and B. Batlogg, Nature 406, 702 
(2000). 

''2 J. H. Schon, C. Kloc, and B. Batlogg, Nature 408, 549 
(2000). 

"3 Y. Maeno, T. M. Rice, and M. Sigrist, Physics Today 54, 
42 (2001). 

L. J. Buchholtz and G. Zwicknagl, Phys. Rev. B 23, 5788 
(1981). 

^5 C. R. Hu, Phys. Rev. Lett. 72, 1526 (1994). 

''^ K. Sengupta, I. Zutic, H. J. Kwon, V. M. Yakovenko, and 

S. D. Sarma, Phys. Rev. B 63, 144531 (2000). 

A. J. Leggett, Rev. Mod. Phys. 47, 331 (1975). 

G. M. Banner and P. M. Chaikin, Phys. Rev. Lett. 75, 

4690 (1996). 

''^ In the calculation of Xmn {T, H) the full effect of the pres- 
ence of vortices need to be incorporated. A microscopic 
calculation of Xmn {T, H) taking into account the spatial 
inhomogeneities introduced by the magnetic field and vor- 
tex core states is underway. 

5° Vortex core states are not as normal as the normal state 
of the superconductor, they have a lot of structure even in 
conventional singlet "s- wave" superconductors. Vortex core 
states in singlet "d-wave" superconductors have also a rich 
structure. So it is expected that vortex core states in triplet 
"p-wave" superconductors are even richer in structure due 
to the additional spin degrees of freedom present. The vor- 
tex core structure will be discussed in a later publication. 

51 C. Caroli, P. G. de Gennes, and J. Matricon, Phys. Lett. 
9, 307 (1964). 

52 J. Bardeen, R. Kummel, A. E. Jacobs, and L. Tewordt, 
Phys. Rev. 187, 556 (1969). 

53 L. Kramer and W. Pesch, Z. Phys. 269, 59 (1974). 

5-* S. UUah, A. T. Dorsey, and L. J. Buchholtz, Phys. Rev. B 



17 



42, 9950 (1990). 

F. Gygi and M. Schliiter, Phys. Rev. B 43, 7609 (1991). 
C. A. R. Sa de Melo, Phys. Rev. Lett. 73, 1978 (1994). 
M. Ichioka, N. Hayashi, N. Enomoto, and K. Machida, 
Phys. Rev. B 53, 15316 (1996). 

M. Franz and Z. B. Tcsanovic, Phys. Rev. Lett. 80, 4763 
(1998). 

C. A. R. Sa de Melo, Phys. Rev. B 60, 10423 (1999). 
The two fluids are not strictly independent, because vortex 

core states depend on properties (like the order parameter 
behavior) outside and at the boundaries of the core. 



B. S. Shivaram, T. F. Rosenbaum, and D. G. Hinks, Phys. 
Rev. Lett. 57, 1259 (1986). 

C. H. Choi and J. A. Sauls, Phys. Rev. Lett. 66, 484 (1991). 
K. Mortensen, Y. Tomkiewicz, and K. Bechgaard, Phys. 
Rev. B 25, 3319 (1982). 

I. J. Lee and M. J. Naughton, The Superconducting State 
in Magnettc Fields: Special Topics and New Trends (World 
Scientific, Singapore, 1998), chap. 14, pp. 272-295, (Edited 
by C. A. R. Sa de Melo). 



